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Theory and simulation of photogeneration and transport 
in Si-SiO x superlattice absorbers 
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Abstract Si-SiO^ superlattices are among the candi- 
dates that have been proposed as high band gap ab- 
sorber material in all-Si tandem solar cell devices. Due 
to the large potential barriers for photoexited charge 
carriers, transport in these devices is restricted to quan- 
tum confined superlattice states. As a consequence of 
the finite number of wells, large built-in fields and any 
kind of disorder, the electronic spectrum can deviate 
considerably from the minibands of a regular superlat- 
tice. In this paper, a quantum-kinetic theory based on 
the non-equilibrium Green's function formalism for an 
effective mass Hamiltonian is used to investigate photo- 
generation and transport in such devices for arbitrary 
geometry and operating conditions. By including the 
coupling of electrons to both photons and phonons, the 
theory is able to provide a microscopic picture of indi- 
rect generation, carrier relaxation and inter-well trans- 
port mechanisms beyond the ballistic regime. 



Keywords solar cell 
port • NEGF 



1 Introduction 



superlattice • quantum trans- 



Si-SiO x superlattices have been proposed as candidates 
for the high band gap absorber component in all-Si tan- 
dem solar cells [TJ[5] . In these devices, photocurrent flow 
is enabled via the overlap of states in neighboring Si 
quantum wells separated by ultra-thin oxide layers, i.e. 
unlike in the case of an intermediate band solar cell, 

U. Aeberhard 

IEF-5: Photovoltaik, Forschungszentrum Jiilich, D-52425 Jiilich, 
Germany 

Tel: +49 2461 61 2615 

E-mail: u.acberhard@fz-juelich.de 



the superlattice states contribute to the optical transi- 
tions and at the same time provide transport of photo- 
carriers, which makes it necessary to control both the 
optical and the transport properties of the multilayer 
structure. To this end, a suitable theoretical picture of 
the optoelectronic processes in such type of structures 
is highly desirable. 

There are several peculiar aspects of the device which 
require special consideration in the choice of an ap- 
propriate model. First of all, a microscopic model for 
the electronic structure is indispensable, since the rel- 
evant states are those of an array of strongly coupled 
quantum wells. In a standard approach, these states 
are described with simple Kronig-Penney models for a 
regular, infinitely extended superlattice. The superlat- 
tice dispersion obtained in this way can then be used 
to determine an effective density of states as well as 
the absorption coefficient to be used in macroscopic 
ID solar cell device simulators. However, depending on 
the internal field and the structural disorder, the het- 
crostructure states may deviate considerably from regu- 
lar minibands or can even form Wannier-Stark ladders. 
Furthermore, the charge carrier mobility, which has a 
crucial impact on the charge collection efficiency in so- 
lar cells, depends on the dominant transport regime at 
given operating conditions, which may be described by 
miniband transport, sequential tunneling or Wannier- 
Stark hopping [3] , relying on processes that are not ac- 
cessible to standard macroscopic transport models. 

In this paper, the photovoltaic properties of quan- 
tum well superlattice absorbers are investigated numer- 
ically on the example of a Si-SiO^ multilayer structure 
embedded in the intrinsic region of a p-i-n diode, using 
a multiband effective mass approximation for the elec- 
tronic structure and the non-equilibrium Green's func- 
tion (NEGF) formalism for inelastic quantum trans- 
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port, which permits to treat on equal footing both co- 
herent and incoherent transport as well as phonon-assisted 
optical transitions at arbitrary internal fields and hct- 
erostructure potentials. 



2 Theoretical Model 

In order to enable a sound theoretical description of the 
pivotal photovoltaic processes in semiconductor nanos- 
tructures, i.e. charge carrier generation, recombination 
and collection, both optical transitions and inelastic 
quantum transport are to be treated on equal foot- 
ing within a consistent microscopic model. To this end, 
a theoretical framework based on the NEGF formal- 
ism was developed [H[S] and applied to quantum well 
solar cell devices. Here, we reformulate the theory for 
a multiband effective mass Hamiltonian, similar to [51 
[7], and extend it to cover the phonon-assisted indirect 
transitions that dominate the photovoltaic processes in 
Si-based devices. Furthermore, in difference to the for- 
mer case, both photogeneration and transport processes 
take place within superlattice states, since escape of car- 
riers to continuum states is not possible due to the large 
band offsets. 



2.1 Hamiltonian and basis 

The full quantum photovoltaic device is described in 
terms of the model Hamiltonian 



H =H e + 
H e =H° + HI, 



H 



v 



H e H e *y -\- H e p ~\- H ee , 



(total) (1) 
(electronic) (2) 
(interaction) (3) 



consisting of the coupled systems of electrons (H e ), pho- 
tons (-ff 7 ) and phonons (H p ). Since the focus is on the 
electronic device characteristics, only H e is considered 
here, however including all of the terms corresponding 
to coupling to the bosonic systems. 

The electronic system without coupling to the bosonic 
degrees of freedom is described by 



H° e = 
with 



2m 



A+U(z), 



U(z) = U(z) + V {z), 



(4) 



(5) 



where Vo is the heterostructure potential and U is the 
Hartree term of the Coulomb interaction correspond- 
ing to the solution of Poisson's equation that considers 
carrier-carrier interactions (H ee ) on a mean-field level. 



The Hamiltonian representations for the interaction 
terms are obtained starting from the single particle in- 
teraction potentials. For the electron-photon interac- 
tion, the latter is given via the linear coupling to the 
vector potential operator of the electromagnetic field 
A, 



H r 



m 



-Ap 



(6) 



with p the momentum operator and 
A(r,i) =J2 [A (A,q)S A)q (t)+AS(A,-q)St _ q (t) 



A (A,q) = e Aq , 



(7) 
(8) 



where e Aq is the polarization of the photon with wave 
vector q and energy Sw q added to or removed from 
photon mode (A, q) by the bosonic creation and anni- 
hilation operators 



5* 



(9) 



and V is the absorbing volume. 

The vibrational degrees of freedom of the system are 
described in terms of the coupling of the force field of 
the electron-ion potential V e i to the quantized field VI 
of the ionic displacement [5] , 



H ep (r, t) = U(L + k, t) ■ W d [r - (L + «)], (10) 
with the displacement field given by the Fourier expan- 



Z4(Lk, t) =J2U* K (A, Q)e i;Q '( L+K > [a AtQ (t) + a^_ Q (t)] , 

A.Q 

(11) 

where the ion equilibrium position is L + k, with L the 
lattice position and k the relative position of a specific 
basis atom at this lattice site, and fl/i.Q,a^Q are the 
bosonic creation and annihilation operators for a (bulk) 
phonon mode with polarization A and wave vector Q in 
the first Brillouin zone. The corresponding polarization 
vector is e QK7 i(Q). 

For numerical implementation of the model, the above 
Hamiltonian needs to be represented in a suitable basis. 
Due to the amorphous nature of the SiO^ layers, atom- 
istic models are of limited applicability. Furthermore, 
the use of an effective mass theory simplifies the elec- 
tronic model considerably. For a quasi-onedimensional 
multilayer system, where quantization appears only in 
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the vertical (growth) direction, the corresponding basis 
functions have the form 



V'mk,, (r) = (p lk Jr)u nko (r), 



(12) 



where tp%k,, is the envelope basis function for discrete 
spatial (layer) index i (longitudinal) and transverse mo- 
mentum ky , u„k is the Bloch function of bulk band 
n, centered on ko. In the case of a system with large 
transverse extension, the envelope basis function can be 
written as 



tkiirii 



Pik,, (r) 



(13) 



where r|| = (x,y), S is the cross sectional area and 
Xi is the localized longitudinal envelope function basis 
element. For the latter, finite element shape functions 
are a popular choice [TJIH] • Here, we will use a simple fi- 
nite difference basis equivalent to a separate single band 
tight-binding approach for each band [TOTTTIH^] . In the 
above basis, the fermion field operators for the charge 
carriers are represented via 



2,n,k|| 



(14) 
(15) 



j,n,ki 



where , c are single fermion creation and annihila- 
tion operators. The representation of the model system 
Hamiltonian in the above basis is now obtained in stan- 
dard second quantization, i.e. 



H(t) = / <Pr&(r,t)H°$(r,t) 



H *njjfi»(k||)4ik|| (*)Cjmk|| (*). 



(16) 

(17) 



J ,,j n,m ky 



2.2 Green's functions, self energies and quantum 
kinetic equations 

Within the non-equilibrium Green's function theory of 
quantum optics and transport in excited semiconduc- 
tor nanostructures, physical quantities are expressed in 
terms of quantum statistical ensemble averages of sin- 
gle particle operators for the interacting quasiparticles 
introduced above, namely the fermion field operator 
& for the charge carriers, the quantized photon field 
vector potential A for the photons and the ionic dis- 
placement field U for the phonons. The corresponding 



Green's functions are 



^a,2)=-^A(i)i fc (2)) c , 

Vlp(l,2) = -Uu a {l)Up(2))c, 



(electrons) (18) 
(photons) (19) 
(phonons) (20) 



where (■■■}c denotes the contour ordered operator av- 
erage peculiar to non-equilibrium quantum statistical 
mechanics [TBirPf] for arguments 1 = (ri,ii) with tem- 
poral components on the Keldysh contour 14J. 

The Green's functions follow as the solutions to cor- 
responding Dyson's equations |15}I16 |J17"U 8]. 



/ d3 [So" 1 (1,3) - 27(1,3)] 0(3,2) = 5(1-2), 



J da^r^a) -^a.3)] w t (3,2) = *a-2) 

J d3 [(D p Q )-\l, 3) - WO, 3)] 2> p (3, 2) = 5(1 - 2). 

(21) 

Go, an( i are the propagators for noninteract- 
ing electrons, photons and phonons, respectively, O de- 
notes transverse and boldface tensorial quantities. The 
electronic self-energy S encodes the renormalization of 
the charge carrier Green's functions due to the interac- 
tions with photons and phonons, i.e. generation, recom- 
bination and relaxation processes. Charge injection and 
absorption at contacts is considered via an additional 
boundary self-energy term reflecting the openness of 
the system. The photon and phonon self-energy tensors 
^T* 7 and II P describe the renormalization of the optical 
and vibrational modes, leading to phenomena such as 
photon recycling or the phonon bottleneck responsible 
for hot carrier effects. The self-energies can be derived 
cither via perturbative methods using a diagrammatic 
approach or a Wick factorization or using variational 
derivatives. In the following, any renormalizing effect 
of the electronic system on the photons and phonons is 
neglected, i.e. the coupling to the bosons corresponds to 
the connection to corresponding equilibrium reservoirs. 
While this treatment is generally a good approximation 
in the case of phonons, it is valid for the coupling to the 
photonic systems only in the case of low absorption, i.e. 
weak coupling or very short absorber length. 

The use of the equilibrium boson propagators im- 
plies that only the electronic Dyson equations are solved. 
In the chosen discrete real-space basis, the components 
of the steady-state Dyson and Keldysh equations for 
electronic Greens functions are turned into a linear sys- 
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ten£J> = ky,E) 

gV) = [{g^)}- 1 -s^h-s^(.) 

G*(v)= [(E + i^l-h^)]- 1 , 
G A (p) = [G»]t, 

G$(v) = G fl (z/) + S^ B 0)1 G A (^) 



(22) 

(23) 
(24) 

(25) 



for each total energy E and transverse momentum k|| . 
There are two types of self-energies in the above equa- 
tions. The term IJ' 1 is due to the interactions between 
the different degrees of freedom of the system. The ex- 
pressions for electron-photon and electron-phonon in- 
teraction are determined as the Fock term within many- 
body perturbation theory on the level of a selfconsistent 
Born approximation, and using the equilibrium boson 
propagators, they are obtained in the following form, 
(a = 7,p) 

Eg (ky ; E) = M ea (k,| , q, A) [N^G$ (k||5 E T ^A, q ) 

A,q 

+ (iV A Q q + l)G^(k||;£;±^ q )] 

x M ea (kii,-q,A) (26) 



and 



£f/(k,|;£)=i 



V 



dE' 17> (,k||;gQ-17< a (k||;^) 
2ir E' -E±ir) 

dE' r(\s. {l ;E') i 
E'-E T 2 



T-r(k f E), (27) 



where 



r(k|, ;E)=i [E> a (k,| ; E) - 17< (ky ; E)] 



(28) 



Since the principal value integral corresponds to the 
real part of the self energy and thus to the renormal- 
ization of the electronic structure, which is both small 
and irrelevant for the photovoltaic performance, it is 
neglected in the numerical implementation. 

Once the Green's functions and self-energies have 



been determined via self-consistent solution of Eqs. ( 22 ) 



(25) and (26)- (27 1, they can be directly used to express 



the physical quantities that characterize the system, 
such as charge carrier and current densities as well as 
the rates for the different scattering processes. 

1 In steady-state, the Green's functions depend only on the 
difference r = t — t' of the real-time variables, which is Fourier- 
transformed to energy. 



2.3 Microscopic optoelectronic conservation laws and 
scattering rates 

The macroscopic balance equation for a photovoltaic 
system is the steady state continuity equation for the 
charge carrier density 



V-j c (r) =Q c (r) -TZ c (r), c = e,h, 



(29) 



where j c is the particle current density, Q c the gener- 
ation rate and 1Z C the recombination rate of carriers 
species <Q In the microscopic theory, the divergence of 
the electron (particle) current is given by [TBirH] 



V-j(r) = - 



dV 



E R (r,r';E)G<(r\r;E) 



dE 
V J "hxh _ 
+£< (r, r'; E)G A (r / , r; E) - G R {r, r'; E) 

xZ<(r',r;E)-G<(r,r';E)£ A (r',r:E) 

(30) 

If the integration is restricted to either conduction or 
valence bands, the above equation corresponds to the 



microscopic version of (29) and provides on the RHS 



the total local interband scattering rate. The total in- 
terband current is found by integrating the divergence 
over the absorbing/emitting volume, and is equivalent 
to the total global transition rate and, via the Gauss 
theorem, to the difference of the interband currents at 
the boundaries of the interacting region. Making use of 
the cyclic property of the trace, it can be expressed in 
the form 



d A r 



dE 
2ttH 



d\' 



S < (r,r';E)G > (r',r;E) 



S > (r 7 r';E)G < (r',r;E) , 



(31) 



_1 . If we are interested in the in- 



with units [R] = 
terband scattering rate, we can neglect in Eq. (31) the 



contributions to the self-energy from intraband scat- 
tering, e.g. via interaction with phonons, low energy 
photons (free carrier absorption) or ionized impurities, 
since they cancel upon energy integration over the band. 
Since inequivalent conduction band valleys may be de- 
scribed by different bands, the corresponding inter-valley 
scattering process has also interband character with a 
nonvanishing rate, as long as only one of the valleys 
is considered in the rate evaluation. Furthermore, if 
self-energies and Green's functions are determined self- 
consistently as they must in order to guarantee current 
conservation, the Green's functions are related to the 
scattering self-energies via the Dyson equation for the 
propagator and the Keldysh equation for the correlation 



functions as given in Eqs. (22)-(25), and will thus be 



The dimensions are that of a volume rate, [5,7?-] 
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modified due to the intraband scattering. In the present 
case of indirect optical transitions, the Greens func- 
tions entering the rate for electron-photon scattering 
between the r bands are the solutions of Dyson equa- 
tions with an intcrvalley phonon scattering self-energy 
and may thus contain contributions from the X-valleys. 
In the same way, the r c Greens functions entering the 
electron-phonon r c — X scattering rate contain a photo- 
generated contribution. By this way, indirect, phonon- 
assisted optical transitions are enabled. 

3 Implementation for Si-SiO x superlattice 
absorbers 

3.1 Electronic structure model 

Within the EMA for silicon chosen for this work, the 
electrons are described by a multi-valley picture with 
different values for transverse and longitudinal effec- 
tive mass, similar to |18j . However, for simplicity, in 
the case of transverse X valleys (X\\), the anisotropy in 
the transverse mass is neglected and an average value 
is used. The virtual r states used in the indirect tran- 
sitions are described by an additional (negative) mass. 
The holes are modelled by two decoupled single bands 
with different effective masses corresponding to heavy 
and light holes. Thus, in total, five bands are used 
to describe the electronic structure, three for the elec- 
trons (Xii, X±, r c ) and two for the holes (-T„;, r v h). 
The band parameters used in the simulations are listed 
in Tab. [T] For each band, a set of Green's functions 
are computed from the corresponding decoupled Dyson 
and Kcldysh equations. In the computation of physical 
quantities such as electron and hole densities as well 
as the corresponding current densities, the summation 
over all conduction or valence bands needs to be per- 
formed, 

rii = ^ f*> n i,b (32) 
b=r c ,x a ,x ± 

= / ^(-^ )6 (k||;£), (33) 

b k|l 

where ft, denotes the degeneracy of the conduction bands, 
which is fr c = 1, fx« = 4 and fx ± = 2. Similarly, the 
electron current in terms of the Green's functions reads 

^ = ]T f b Ji, b (34) 
b =r a ,x ll ,x ± 

b kg J 

-k + iiGf i+ljb (kfE)]. (35) 



For the chosen model of the bulk band structure, the 
total radiative rate is 

k n 

-£<. rc (k||;£)G>>|| ;£)]}, (36) 
and the inter-valley phonon scattering rate reads 

R ^r-x = ~l f Tr{£ [*> p(r _ x) , rc (h;E) 

k|, 

x G< (k„ ; E) - S< {r _ x) rc (k„ ; E)G> e (k„ ;£)]}. 

(37) 



3.2 Interaction parameters 

Optical transitions are assumed to take place only at 
the center of the Brillouin zone, i.e. between r v and 
virtual F c states, the latter being (de-)populated via 
phonon scattering from (to) the X valleys, which carry 
the photocurrent. All other transition channels, e.g. phonon 
scattering in the valence band prior to photon absorp- 
tion, are neglected at this stage. The momentum ma- 
trix element in the electron-photon coupling is thus to 
be taken between the r v and r c bands at ko = 0. The 
interaction matrix elements are evaluated using an av- 
erage effective coupling for both light and heavy holes. 

Four different types of phonons are used in the present 
work to describe both carrier relaxation as well as phonon 
assisted optical transitions. For the relaxation process, 
X — X inter-valley scattering is used for the electrons 
and non-polar optical phonon scattering for the holes. 
Further broadening is added for both carrier species 
through acoustic phonon scattering in the deformation 
potential formulation. Finally, the momentum trans- 
fer for the indirect optical transitions is mediated via 
r c — X intervalley scattering. 

Table 1 Band parameters used in simulations 





Si 


SiO* 


m* rc /mo 


-0.3 


-0.3 


m *x n / m o 


0.98 


0.4 




0.19 


0.4 


m *rv,ih/ m o 


0.16 


0.4 


m rv,hJ m o 


0.49 


0.4 


E g: rv-ra [eV] 


3.5 


5.5 


E g ,rv-x [eV] 


1.1 


3.1 
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intrinsic absorber 




Fig. 1 Spatial structure of the p-i(SL)-n model system. 

4 Numerical results and discussion 

4.1 Model system 

The model system under investigation is shown schemat- 
ically in Fig. [I] It consists of a set of four coupled quan- 
tum wells of 6 monolayer (ML) width with layers sep- 
arated by oxide barriers of 3 ML thickness, embedded 
in the intrinsic region of a Si p-i-n diode. The thickness 
of the doped layers is 50 ML, while the total length of 
the i-region amounts to 154 ML. This composition and 
doping leads to the band diagram shown in Fig. [2] 

4.2 Density of states 

Insertion of the oxide barriers leads to an increase of 
the effective band gap in the central region of the diode 
from 1.1 eV to ~ 1.3 eV, as seen in Fig. [3j which shows 




1 1 1 1 11 1 ^ 

10 20 30 40 50 60 



position z [nm] 

Fig. 3 Transverse momentum integrated local density of states 
of the p-i(SL)-n photodiode at short circuit conditions. 




position 



Fig. 2 Band diagram of the p-i(SL)-n model system. 

the transverse momentum integrated local density of 
states. In the actual situation of strong band bending, 
quantization also occurs in the form of notch states in 
front of the barriers. The density of states at minor- 
ity carrier contacts is additionally depleted due to the 
imposition of closed system boundary conditions that 
prevents the formation of a dark leakage current under 
bias. The density of states component at zero transverse 
momentum displayed in Fig. [4] allows the identification 
of the confined states in the different quantum wells, 
which are considerably localized due to the large inter- 
nal field, with however finite overlap between neighbor- 
ing wells in the case of the higher states. The ground 
state is split due to the different effective masses of the 
charge carriers, the effect being more pronounced for 
the electrons. 




1 1 LJ - t 1 1 ' !— 1 1 

20 25 30 35 40 45 

position z [nm] 

Fig. 4 Local density of states in the quantum well region at zero 
transverse momentum (kii = 0). 
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Fig. 5 Spatially and energy resolved photon carrier photo- 
generation rate in the quantum well region at short circuit 
conditions and under monochromatic illumination with en- 
ergy E-y = 1.65 eV and intensity / 7 = 10 kW/m 2 . 

4.3 Generation and photocurrent spectrum 

The spectral rate of carrier generation in the confined 
states under illumination with monochromatic light at 
photon energy E 1 = 1.65 eV and intensity 7 7 = 10 
kW/m 2 is shown in Fig. [5] At this photon energy, both 
lowest and second minibands are populated. The pho- 
tocurrent originating in this excitation is shown in Fig. 
[H Current flows also in both first and second mini- 
bands, i.e. over the whole spectral range of generation, 
which means that relaxation due to scattering is not 
fast enough to confine transport to the band edge. How- 
ever, transport of photocarriers is strongly affected by 
the inelastic interactions and is closest to the sequential 
tunneling regime. 

5 Conclusions 

In this paper, an adequate theoretical description of 
photogeneration and transport in Si-SiO^ superlattice 
absorbers was presented. Based on quantum kinetic the- 
ory, the formalism allows a unified approach to both 
quantum optics and inelastic quantum transport and 
is thus able to capture pivotal features of photogener- 
ation and photocarrier extraction in Si-based coupled 
quantum well structures, such as phonon-assisted op- 
tical transitions and field-dependent transport in su- 
perlattice states. Due to the microscopic nature of the 
theory, energy resolved information can be obtained, 
such as the spectra for photogeneration rate and pho- 
tocurrent density, which shows that in the case of high 
internal fields, excess charge is transported via sequen- 
tial tunneling in the miniband where it is generated. 




28 30 32 34 36 

position z [nm] 



Fig. 6 Spatially and energy resolved photon carrier short- 
circuit photocurrent density in the quantum well region under 
monochromatic illumination with energy E 7 = 1.65 eV and 
intensity / 7 = 10 kW/m 2 . 
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